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Abstract 

Stationary and axially symmetric space-times play an important role in astro¬ 
physics, particularly in the theory of neutron stars and black holes. The static 
vacuum sub-class of these space-times is known as Weyl’s class, and contains 
the Schwarzschild space-time as its most prominent example. This paper is 
going to study the space of Killing tensor fields of valence 3 for space-times 
of Weyl’s class. Killing tensor fields play a crucial role in physics since they 
are in correspondence to invariants of the geodesic motion (i.e. constants of the 
motion). 

It will be proven that in static and axially symmetric vacuum space-times 
the space of Killing tensor fields of valence 3 is generated by Killing vector fields 
and quadratic Killing tensor fields. Using this result, it will be proven that for 
the family of Zipoy-Voorhees metrics, valence-3 Killing tensor fields are always 
generated by Killing vector fields and the metric. 

1 Introduction 

Consider a manifold M with Lorentzian metric g , and its cotangent space T*M 
endowed with its natural symplectic form. A Killing tensor field K of valence d 
on M is a symmetric (0, d)-tensor such that 

V(aK bl ...b d ) = 0 - ( 1 ) 

The Lorentzian metric provides an isomorphism between T*M and TM, and 
we are therefore going to identify co- and contravariant tensor fields as well as 
the corresponding homomorphisms with mixed co- and contravariant indices. 

Killing tensor fields are in 1-to-l correspondence to (first) integrals (or 
Hamiltonian invariants) polynomial in the momenta p (or the velocities 7 ) of 
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the Hamiltonian motion for the Hamiltonian function H = g ZJ PiPj = g( 7 , 7 ). In 
the language of integrals, the Killing tensor equation takes the form 


{Ik,H} = X h {Ik ) = 0 


( 2 ) 


where {•, •} denotes the usual Poisson bracket on T*M, and where Xh(Ik) is the 
derivative of the function Ik = K( 7 ,... , 7 ) in the direction of the Hamiltonian 
vector field X h- Two integrals are said to be in involution if they commute 
with respect to the Poisson bracket 

Studying the existence of polynomial integrals is interesting from at least two 
perspectives. Firstly, the existence of integrals can help in answering natural 
questions about the behavior of trajectories, i.e. the behavior of free falling parti¬ 
cles in physically motivated Hamiltonian systems (e.g. by the famous Liouville- 
Arnold Theorem one can integrate the system by quadratures under certain 
additional assumptions). 

Secondly, asking for the existence of integrals is a natural geometric re¬ 
quirement. Metrics meeting this requirement may lead to physically interesting 
examples, as for e xampl e in the case of the Kerr metric possessing the Carter 


constant 


Car 68 a : Car 68 b| | . an integral quadratic in momenta in addition to 
energy and axial symmetry. Integrals polynomial in the momenta are of par¬ 
ticular interest since they represent a generalization of constants of the motion 
that emerg e from the action of one-parameter groups of diffeomorphisms. Ac¬ 
cording to Mar 14], the example of Kerr-de Sitter space-times is at present the 
only known example of integrable space-times with an additional integral of 
higher-than-linear degree, among the class of stationary and axially symmetric 
space-times. 

Static and axially symmetric vacuum (StAV) space-times form Weyl’s class ; 
they a re spec ial cases of stationary and axially symmetric vacuum (SAV) space- 
times ISte03l : BriOSb]. Recently, some attention has been drawn to the Zipoy- 


Voorhees family , which belongs t o this class Voo7 0: Zip 66 ]. Numerical studies 
Bri08a : Bri08b; BrilOat iBril Ob] suggested integrabili ty for some StAV met - 


rics, while later studies provided contradicting evidence KM12 : Lukl 2 t MPS13 |. 
Note that while these studies considered fixed values of the parameter <5 of the 
Zipoy-Voorhees metric, in this paper we are going to consider arbitrary 6 in the 
case of the Zipoy-Voorhees metric. Therefore the result for the Zipoy-Voorhees 
family is in line with the evidence contradicting integrability of the family. 

The methods used in this paper are not restricted to the concrete setting 
of Weyl’s class. It is therefore likely that the methods are suitable for other 
examples of parametrized metrics in two dimensions. It seems probable that 
such examples include 2 -dimensional ones with potential, since this is closest to 
the case studied here. The study of integrals in 2-dimensional manifolds is a 
classical problem in differential geometry and goes at least back to Darboux and 
Koenigs [Dar87i] . For instance, some non-existe nce res ults have been obtained 
on the 2-torus for cubic and quartic integrals in Bva87t NVDOOl : BM11 1, while 
for higher degrees almost nothing is known. On the 2-sphere, however, some 
integrable examples are known. For instance, a new integrable system has been 
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presented, with the additional integral being of cubic degree, by Dullin and 
Matveev DM041 ]. 

The procedure taken in this paper is a new approach to the question of 
existence of integrals. It combines two major previous lines of action: 

- It is inspired by ideas from prolongation-projection methods well-known 
in the theory of overdetermined PD E syste ms. However, it does not follow 
the algorithmic procedure used in KM12 |. 


- It follows ideas outlined in Hie87j] . but takes a somehow converse track 
that could be described as “bottom-up” in contrast to the “top-down” 


approach taken in Hie87l |. The advantage of this direction of reasoning 
is that it avoids solving the leading-degree Poisson equation. Instead, 
our approach rather begins with solving simple geometric orthogonality 
relations. 


1.1 Static and Axially Symmetric Vacuum Space-Times 

Stationary and axially symmetric vacuum space-times possess two commuting 
Killing vector fields, one being space-like and the other being ti me-lik e. Such 
space-times can be brought into the following standard form SteOS ] by the 
aid of suitable coordinat e trans for matio ns. The coordinates are called Lewis- 
Papapetrou coordinates fBri08a ; BrilObl ]. 


g = e 2U (e 27 (dx 2 + dy 2 ) + R 2 d</> 2 ) — e 2U (dt — oj d(f>) 2 


(3) 


We will restrict our attention to vacuum space-times and therefore require the 
Ricci tensor of g to be identically zero. This is a fair assumption for the move¬ 
ment of test particles around astrophysical objects as long as electro-magnetic 
fields are ignored. For SAV space-times, Ricci flatness is encoded in a set of 
equations which are called the Ernst equations. In the static case, we require 
to = 0. Then, the Ernst equations read as follows 


Ry Uy + R x U x + R Uyy + R U xx — 0 
AR = R xx + Ryy = 0 

2 R U 2 — 2 R y 7 y + 2 R x ") x + R xx — R yy — 2 R U 2 = 0 
2 R U x U y + R y 7 X + R x 7 y + R xy = 0 

The equations break up into two sets of two equations each, which we shall 
refer to as primary and secondary equations. The primary equations give re¬ 
strictions on U and R. Provided R is non-c onstant, A R = 0 allows setting 
R = x > 0 by a change of coordinates [Ste03l ]. If R is constant, this change of 
coordinates is impossible, but one can show that A 7 = 0 holds and the metric 
is flat. In case of non-constant R , the secondary equations enable us to express 
derivatives of 7 in terms of derivatives of U , allowing us to eliminate them, and 
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finally 7 , from the equations. We obtain the relations: 


Uyy TJXX U x 

R = x 


7 x = -x U 2 +xU 2 

'fy — ^.x U x Uy 


(4) 


By definition, stationarity and axial symmetry can be described by the global 
symmetry group h = 1R, x S 1 . The action of h is Hamiltonian and we denote the 
moment map by p. h acts freely, and we may pass to the symplectic quotient 
Qred = /U _1 (0)//i, which inherits a symplectic form from the initial space-time 
(with an additional com pactness assumption, this is the Marsden- Weinstein 
quotient , see e.g. [MS95 |). In Lewis-Papapetrou coordinates, h acts along coor¬ 
dinate directions and we will be able to identify the reduced coordinates easily. 
The 4-dimensional problem then is reformulated as a 2-dimensional problem 
with metric <?red> and the Hamiltonian H = T + V splits into a kinetic term 
T = H led = g\ J ed PiPj along with a potential V, which is polynomial in p$ and p t . 
As a consequence, we distinguish Hamiltonian integrals that commute with H, 
from metric integrals that commute with T only. Note that the highest-degree 
component w.r.t. ( p x ,Py ) of a Hamiltonian integral is metric (i.e. commutes with 
T). The metric and the potential on the reduced space read: 


5 red = e 2U (e 27 (dx 2 
->- 2 - 2 U „2 


dy 2 )) 


V = R- z e~ 2IJ p 2 -e 2U p 2 


(5a) 

(5b) 


1.2 Main Results 

We consider integrals of the general form 

3 i 3—i 

i (x,y) = Yl a i,xk( x ,y) pip i y~ j p\p\~ r ~ k ( 6 ) 

2—0 j—0 k—0 

Such integrals are in involution with p ^ and pt . Since the Hamiltonian defined 
by equations m does not mix momenta p x ,Py with p$,pt, the components 
of © of odd and even parity in (p x ,p y ) can be considered separately. We prove 
reducibility of degree-3 integrals in space-times of Weyl’s class. 

Definition 1. Let I be a polynomial integral of degree d. We say that I is 
reducible (by one degree) if there are polynomial integrals Ji,... ,/ m of degree 
at most d — 1 such that I is a linear combination of products of the integrals R. 
We say that I is totally reducible if there is a representation of this form such 
that the Ii are integrals linear in momenta or the Hamiltonian. 

For space-times in Weyl’s class, we prove reducibility of degree 3 integrals 
by one degree. In addition, total reducibility of degree-3 integrals is shown for 
the family of Zipoy-Voorhees space-times (a sub-family of Weyl’s class). 

Theorem 1. Let M be a space-time in Weyl’s class. Then any integral © of 
third degree on M is reducible. 
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Reducibility of a degree-3 polynomial integral means that these integrals 
can be written via products of lower-degree integrals. In the language of Killing 
tensors, this means that any valence-3 Killing tensor field can be written via 
symmetrized products of Killing vector fields and quadratic Killing fields. 

For a concrete example, consider the Zipoy-Voorhees class of space-times. 
Their metrics are static and axially symmetric, and parametrized by a parameter 
S > 0 Voo70l |. Therefore, this family forms a subset of Weyl’s class. 



+ (X 2 -1)(1 (7) 

The resulting metric for S = 0 is flat. The value <5 = 1 gives the Schwarzschild 
metric (that admits one additional quadratic integral). We allow arbitrary S > 
0 . 

Proposition 1. Let Mzv be a Zipoy-Voorhees metric with parameter 6 > 0, 
5^1. Let I be an integral © of third degree on Mzv- Then I is totally 
reducible, i.e. generated by linear integrals (i.e. Killing vector fields) and the 
Hamiltonian (i.e. the metric). 


2 Method 

The basic procedure is as follows: 

(i) Reduce the 4-dimensional problem without potential to finding integrals 
on a 2-dimensional Hamiltonian manifold with potential (symplectic re¬ 
duction). 

(ii) The existence of integrals is encoded in equations that emerge from the 
Poisson equation {H, Ik} = 0 as coefficients w.r.t. momenta. Splitting 
according to the degree in momenta ( p x ,Py ) yields three polynomials in 

and pt- If we decompose these polynomials further w.r.t. momenta, we 
obtain three blocks of equations. 

(iii) Use the equations obtained from zeroth degree in momenta (p x ,Py) and 
solve them as far as possible, obtaining one function a to parametrize the 
integral (this is the case without an additional integral present. If there 
is an additional linear integral, lemma [3] applies). 

(iv) From the block obtained from the degree-2 polynomial, extract two inte- 
grability conditions for a. 

(v) The remaining system of equations is an overdetermined system of PDE 
involving the metric which is described by one function U in two variables. 
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We consider derivatives of U as being new, independent unknowns. By 
taking derivatives (prolongation), and then eliminating higher derivatives 
of U (projection), we end up with an ordinary differential equation. 

(vi) For the remaining ODE we show that its only solution corresponds to flat 
space. This allows to conclude that degree-3 integrals are always reducible. 


Equation © is the condition for a function I to be an integral. Since we 
take / to be polynomial in momenta of degree d, © is a polynomial in momenta 
of degree d + 1. We are going to consider the system of PDE obtained from 
the coefficients of © w.r.t. momenta. Symplectic reduction w.r.t. the symmetry 
group (stationarity, axial symmetry) suggests to regard p^ and pt as parameters. 
We distinguish the equations of the PDE system by the momenta monomials to 
which they appeared as a coeffic ient. For Weyl’s class, the equations can then 
be arranged in a tree-like Bri08b| structure. We write down the Hamiltonian in 
the following form: 

H = T+ V^pl + V u p 2 (8) 


=v 


where T = H re ^ is the reduced Hamiltonian (i.e. a homogeneous polynomial 
in p Xl p y , and where V ab are the smooth coefficient functions of p a Pb (a, & £ 
{(/>, t})). The integral I can be decomposed accordingly. We denote 


/=/ <d> +/r > w+4 <d - i, p.+/s 2, f 


T fd-2) 




=/(d-2) 


where each I is of degree k in the momenta p x ,p y . We require the metric to 
be non-flat such that we can choose coordinates with R = x. In this case we 
have three blocks of equations coming from the respective polynomials (this is 
step (ii) of the above list). We can extract the equations from the polynomials 
which are obtained by splitting © according to the degree w.r.t. ( p x ,P y )'■ 


M 

v -v-' 

II 

o 

degree 4 

(10a) 

{T,jW} + {F,/ (3) } = 0 

degree 2 

(10b) 

(E,/ (1 )} = 0 

degree 0 

(10c) 


The equations of even parity in (p x ,p y ) split off from this system and form a 
separate, decoupled system. Equation (llOal) is the condition that must hold for 
an integral I^ on the reduced manifold with Hamiltonian T = iF re d- However, 
only some of these integrals ascent to integrals upstairs on the initial manifold. 
This is due to the restrictions (llObl) and (IIPel) . For a better understanding 
of these equations (or the equations obtained from them as coefficients w.r.t. 
momenta), we will characterize them as defining equations for I^ and /W. 
But first let us split the system further by considering coefficients w.r.t. ( Pt,Pc />)• 
The polynomial (II Hall does not split since it is already of degree 4 in momenta 
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(p x ,Py)- The polynomial (llObl) splits into three parts: 

{t,i$} + {v++,iW} = 0 

0 = {T,I$} 

{T,4P} + {V u ,lW} = 0 

We write the equations in this form to hint at the fact that the equations can 
be divided into two groups that can be treated separately. This will become 
clear when we include (fioa . The second of the equations says that is a 
metric integral on the reduced space. In fact, we will see that it even has to 
be an integral on the initial space-time, and therefore is not of interest for our 
considerations. The polynomial (llQcl) splits into five parts: 

{^4?} = ° 

o = 0^,0 

{V tt ,I$} + {V++,iW} = 0 

o = 

{V U ,4P} = 0 


The second and fourth of these relations tell us that is an integral not only 
on the reduced, but also on the initial space. We can isolate this subsystem 
from the remaining one and solve it separately (this procedure is possible in 
general for Weyl’s class). This can easily be done and is equivalent to finding 
Killing vector fields of the space-time under consideration. 

The remaining equations from w can be interpreted in a nice way as 
scalar product relations for the components of /W. For instance, 

{V H , 1^} = b f 0 + V^bi^ = e 2U ~ 21 (W 00 , b H ) = (,dVb 

where 1^ = bf^p x + b^Py and where \/V^ denotes the gradient vector corre¬ 
sponding to the differential dV^. The polynomial (II Pel) therefore gives rise to 
a set of scalar product relations: 

(VF^, b^) = 0 
(VF M , b H ) + (W w , b u ) = 0 
(VV tt ,b tt )= 0 


This allows to solve (llQcl) directly for b^ and b u , 

b^ = ai V-'-VW, b u = a 2 V-'-V" 

where we introduce the shorthand notation V x / = e 2U ~ 2l (—f y , f x ) for a func¬ 
tion /, i.e. V x / is the vector field rotated by 11 / 2 compared to V/. Defining the 
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angle between VF M and S/V^, 

(S7V U 

COS ^ = --- 

IIVV^II ||W tt ||' 

the second of the three scalar product relations can be brought into the form 

(«2 — Qii) sin'h = 0 (11) 

This is step (iii) in the list given at the beginning of this section. We summarize: 

Lemma 1. Either the metric potentials are such that W^ and W M are 
parallel, or the parameter functions u\ and 0.2 are equal. 

We now turn to an interpretation of (IIObi) . Consider {V,/^} and denote 
j(3) _ pjkp. p . pk Then, 

{F,/ (3) } = 3 C V x r ij p iPj + V y P ij p iPj ) 

= 3(V k I kij PiPj ), 

and analogously for {V^, 1^} and {V u , I^}. With this in mind, we interpret 
(llObl) as defining equations for the tensor field A'( 3 ) (VV, •, •). There are two 
more equations than components of A^ 3 ) and this allows us to find independent 
expressions for A^VV^, •, •) as well as K( 3 \VV U , ■, ■). Then, if dV tt and 
dV^ arc linearly independent, we can determine A'^ 3 ^ in terms of derivatives 
of the function a = on = a^- We have the following obvious identities: 

A (3) (VV^; VF^, W tt ) = A (3) (W m ; WV h> ) (12a) 

A (3) (W^; W M , W“) = A' (3) (W M ; W^ 1 , W M ) (12b) 

We proceed as follows: 

1. Determine A ^ in terms of a and its derivatives, if sin ^ 0. 

2. Determine derivatives of a using the symmetry in the arguments of K^ 3 \ 
Then derive an integrability condition for a. 

3. Combine the integrability condition with (II OaP and the Ernst equations. 
Show that the system does not have any solutions, using algebraic manip¬ 
ulations as well as prolongation-projection arguments for the system. 

First, however, consider the case sin ’I' = 0. In order to do this, we begin with 
a closer look at Killing vectors. 

2.1 Killing Vector Fields 

Assuming there is an additional linear integral on the 4-dimensional space-time, 
we characterize the existence of linear integrals in terms of the rank of the 2x3- 
rnatrix whose columns are given by gradients of the potential components V^, 
V tlt> and V u : 

M = (dV tut> t dV tt ,dV ul '). 
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Since the rank of AA is the dimension of the linear space spanned by dV^, dV tt , 
dV^, it is a geometric object and independent of the choice of coordinates. 

If dV^ = 0, M might be replaced by the 2 x 2-matrix (dV^, dV u ), which 
will also be denoted by AA. Then, instead of the rank of the matrix, the de¬ 
terminant may be used with the obvious correspondences. In case there is an 
additional linear integral present in a non-flat SAV space-time (i.e. in addition 
to and pt), the rank of AA cannot be full. More precisely: 

Lemma 2. (a) Let (M,g) be in the SAV class. 

• If there is an additional linear integral (Killing vector field), then the 
rank of AA is 1, or the space-time is flat. 

• Let rk(A4) = 1. Then p y is a linear integral (Killing vector field) 
when using Lewis-Papapetrou coordinates with R = x. 

(b) Let ( M,g ) be in Weyl’s class. 

• Let rk(A4) < 1 be constant. Then there is an additional linear integral 
on M. In case rk(AI) = 1 this vector field corresponds to p y in 
Lewis-Papapetrou coordinates with R = x; in case rk(_Ad) = 0 the 
space-time is flat. 

• If there is an additional linear integral, it is given by p y in Lewis- 
Papapetrou coordinates with R = x, if M is non-flat. 

Proof. Part (a). For linear integrals we only have two polynomials after taking 
coefficients w.r.t. (p x ,p y ) (they are similar to the polynomials of degree 0 and 2 
in JH 8 ). Let us denote the components of the (p x , p, ; )-linear part of the integral 
as 

/W =b 1 p x +b 2 p y , b=(bi,b 2 ) (13) 

The zeroth order equations read: 

(VV**,b) = 0, (V^,!)) = 0, (VV tt ,b) = 0 (14) 

We conclude that the following relations must hold: 

(vv**,v- L v**) = o, iyv u , = o, (v v**, = o (15) 

This means, the potential gradients are pointing all in the same direction, i.e. 
they are pairwise linearly dependent. Hence, the rank of the potential gradi¬ 
ent matrix is 1, provided the space-time is non-flat. In the flat case we may 
have Lewis-Papapetrou coordinates with the parameter R constant, making it 
impossible to choose R = x. We therefore exclude the flat case from our con¬ 
siderations. This concludes the proof of claim one and establishes a necessary 
criterion for the existence of Killing vector fields in everywhere non-flat space¬ 
times. 

Now, since the rank of the potential gradient matrix AA is 1, all rows as well 
as all columns have to be linearly dependent. This again gives us relations (ED, 
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meaning that V V^, and Vh if are pairwise linearly dependent. First, 

let us assume w/0. We consider 

= o 


This equation amounts to the requirement 

X U! X Uy — (1 + X U X ) U!y = 0 


or the relation 


UJy 


= K 


1 + X U x 

X Uy 


with a scalar function n to be determined. Inserting this into the requirement 

= o 


yields the relation 

U y x 2 e 4U = 0 

and forces U to be a function of x only. Turning back to the relations for u>, we 
see that oj v = 0, so u also is a function of x only. 

Recalling the convention R = x, the metric depends on x only, if 7 only 
depends on x, or if it is constant. We infer the secondary Ernst equations, 

4x 2 e 4U U 2 — u> 2 — 4x e 4U = 0, 

7v x e 4U = 0 . 

Consider the latter equation. It means 7 y = 0, so we are done. Since the metric 
does not depend on y, p y must be an integral, and therefore provides a Killing 
vector field. 

Now, assume u> = 0. Then (W^, S7 ± V t ^ > ) = 0 trivially and we have to go 
a slightly different way of reasoning. We consider 

(WV <M> ,W ± V tt ) = 0. 


It follows that 

x 2 e 4U U y = 0, 

which means U y = 0 (on the entire neighborhood). Conclude U = U(x), and 
then 7 = 7 ( 2 :). Thus, the metric is a function of x only and p y is a linear 
integral. This concludes the proof of part (a). 

For part (b) let us first remark that if an additional linear integral exists in 
Weyl’s class, it must be a multiple of p y , or the metric is flat. Two cases need to 
be checked: Firstly, if there is exactly one additional linear integral, it is a mul¬ 
tiple of p v . Secondly, if there are two (independent) additional linear integrals, 
there are three (say b ^, k = 1, 2,3). Looking at the equations (V V*- 7 ', b^) = 0, 
this forces all gradients to be zero (or, equivalently, dV*- 7 ' = 0). Hence, 
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V is constant. Thus U and ui are constant, in contradiction to the assumption 
R = x. Therefore, the metric is flat. 


With this remark, the first claim of part (b) follows immediately from 
part (a), keeping in mind that rank 0 corresponds to flat space. The second 
claim of part (b) follows immediately from the second statement of part (a). □ 

3 Proof of the Main Theorem 

For the proof we will w.l.o.g. assume constant rank for the matrix A 4. If the 
space-time M is not of constant rk Ai, we may still consider the subsets of points 
in M with constant rank 0, 1, or 2. We may then work with the sets of their 
inner points ignoring the remaining points of M, which amount only to a null 
set w.r.t. the measure induced by the volume form on M. Proving the theorem 
on a dense set is sufficient because if a degree-3 polynomial integral is identical 
to a linear combination of products of H , p$ and pt on an open subset, this is 
true everywhere. 

The proof will be completed in two steps. First we consider space-times with 
rk A4 = 1, then the case rkA4 = 2. As we have seen, the rank-1 case is the 
case when there is one additional Killing vector field. Rank-2 is the case if no 
additional Killing vector field exists (assuming non-flatness). 

Lemma 3. IfrkM. = 1, then any third degree integral is reducible by at least 
one degree. 

Proof. By the hypothesis, there is the linear integral p y in Lewis-Papapetrou 
coordinates with R = x. Consider (II Obi) 

{V++,fW} + {T,F$} = 0 

{V t^ F (3) } + {TjF W }=0 

{V tt ,F^} + {T,F^} = 0 
Each F^ is a multiple of p y , so 

F U = h iPyi F tl ] = h 2Py> F tt ] = h 3 Py 
This means that the equations in (llObj) are of the form 

{V^,F^} + {T, h 1 } Py = 0 
{V t *,FW} + {T,h 2 }p v = 0 
{V tt ,F^} + {T,h 3 } Py = 0 

The leading order term F^ hence is of the form 

F {3) =Px ((■■■)Py)+fp 3 y ='■ Fp y 
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where the leading p x is because the potential gradients (or, equivalently, the dif¬ 
ferentials dV ab ) have only p x components. The final contribution f p y accounts 
for the fact that (II Obi) only specifies components with at least one p x . 

Now consider (llOal) . 


{T J F^} = {T J Fp y } = {T,F}p y = 0. 

This means {T, F} = 0, so F is a quadratic integral on the reduced space. It 
follows that it can be extended to an integral on the initial space-time, because 
of the fact that 


{V^,F(3)} = {V^,Fpy} = {V^,F}py 

and so on, so we have from (llObl) the equations 

{V++,F} + {T,h 1 } = 0 
{V**,F} + {T, h 2 } = 0 
{V u ,F} + {T,h 3 } = 0 

which makes F = F + hi + h 2 p^ pt + pi a quadratic integral on the initial 
space-time (see remark Q] below). Note that F might still be reducible, but can 
be non-reducible as well. □ 

Remark 1. An even-parity quadratic integral I = p^ + iff ptP</> + 

iff p\ satisfies the polynomial equations 

{T,lW} = 0 (16a) 

{V^,I^} + {T,lff} = 0 (16b) 

{V t +,I(*)} + {T,I% ) } = 0 (16c) 

{V tt ,l(*>} + {T,lg ) } = 0 (16d) 

Proof. Decompose {Ft, 1} = 0 by setting each component homogeneous in 
( p Xl p y ) zero. The first equation is the component of degree 3, the other three 
equations are components of degree 1. □ 

We now turn to the case when there is no additional linear integral in invo¬ 
lution with the others. From now on, we will always assume to work in Weyl’s 
class. 

Keeping in mind the considerations of the previous sections, we see that 
this case requires rk.A/1 = 2. Rank 2 requires and V V u to be linearly 

independent. Then, recalling equation m, the scaling functions a± and a 2 are 
equal for Weyl’s class. For simplicity we therefore introduce the new function 
a = a.\ = a 2 (step (iii) in the list of section [2]), so 

b H = aS7 ± V H , b u = a V ± V u . 
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Lemma 4. Derivatives of a are determined by differential equations of the form 

a x = Act 
a v = B a 

where A and B are algebraic expressions determined by Vff, Vff, and Vfj’v, 

which do not contain any higher-than-second derivatives of components of the 
potential V. 

Proof. We use the relations m, i-e- we use the six equations from (II Obi) and 
combine them in a straightforward way to find expressions for the coefficients ao 
through to 0,3 of It = Y^i a iPt~ l Py- I n this way, we find two different expressions 
for ai, and two for 02 , corresponding to the above identities. The expressions are 
polynomials in derivatives of the potential V, i.e. they are determined by Vff. 
Vy f , and and do not contain derivatives of order higher than 2. The 
coefficients of the a,; are just integer multiples of v = (V V u , V 1 - V^), which is 
non-zero because we required V V u and V to be linearly independent. 

We can then eliminate ai and 02 and deduce two equations which have the 
following form: 


(Vy M ,V 1 y^) a x = (...) a 

(vy“,v 1 y # ) a y = (...)a 

The expressions abbreviated by (...) are polynomials in derivatives of V of at 
most second order. Dividing by the non-zero coefficient of the a-derivatives 
yields the required result. □ 

The integrability condition for a is a necessary requirement for the existence 
of non-reducible Killing tensor fields (step (iv) in the list of section [2J : 

Lemma 5. Let rk A4 = 2 and u> = 0, but V V V V tt 7 ^ 0. If there is a Killing 
tensor field of valence 3, then either a = 0 or A y — B x = 0. 

We note that in case a = 0 the integral F 3 = = 0, so the lemma 

actually provides a necessary criterion for the existence of non-trivial Killing 
tensor fields of valence 3. 

Proof of lemma [3 Compute 

(^a:)y (&y)x = A y Ot T A Ot y B x OL B Ot x 

= ( A y — B x ) a + (AB — BA) a 
= i-Ay B x ) ol 

□ 

We give an example where this idea provides information on the reducibility 
of linear integrals: 
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Example 1. The Zipoy-Voorhees family of metrics is a family in Weyl’s class 
that is parametrized by a non-negative number S. 

We can use the method as described above, but we take H in a modified form, 
namely 


H = 


2 fii 


Px 

2fii 


v H v\- 


V tt p 2 t 


The Zipoy- Voorhees metric satisfies, in prolate spheroidal coordinates: 


y<j> 



V* = — 


fx + 1 

— 1 


s 


Taking an approach similar to lemma [ 31 we first check that det M ^ 0. We 
find the following: 


det M = 0 o- 1 5 2 y 2 


x 8 — 4:X 6 y 2 + 6 y 4 a ; 4 — 4 y 6 x 2 + y 8 
{x — l ) 2 (x 2 — l ) 4 (—1 + y 2 ) 6 (x + l ) 2 


= 0 


which obviously is not true for generic x,y, if S ^0. Then we compute the 
necessary criterion as in lemma 0 Using computer algebra (Maple 18), we find 


A y — B x = 0 <=>• 


(x — y ) 4 (x + y ) 4 (—36 x 2 + 4(5 2 x + 2x — 36) 6 2 y i 
(x 2 — l ) 8 (y 2 — l ) 6 


which is not true for generic x,y since 6 > 0. We therefore must conclude 
a = 0, which means the integral of degree 3 is zero. 

We now turn to step (v) in the list at the begining of section [2j 

Lemma 6 . A polynomial equation of degree N > 0 for a function f(x,y) with 
coefficients that depend on x only, is independent of y, so f = f(x). 

Proof. Denote the equation by 'f2n=o cln ( x ) f n ( x ^v) = 0- ^ we can factor out 
f(x , y), then / = 0 and is independent of both x and y. Otherwise, we take one 
derivative w.r.t. y and obtain ^ ra=1 a n (x) n f n ~ x f y = 0. Then either f y = 0 or 
we divide by f y and proceed similarly. At some point, we either get f y = 0 or 
we end up with ajy = 0 , which contradicts the hypothesis that the polynomial 
equation be of degree N. Thus we have f y = 0 and / is a function of x only. □ 

Lemma 7. Let U x = U x (x) be a function of x only. Let the StAV space-time 
have a non-reducible third degree integral. Then U y = 0. 
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Proof. The proof has two parts: (1) Show that U y has to be constant, (2) Show 
that the constant is zero. 

For the first part, consider the p?P 2 -component of (llOal) . Use the Ernst equation 
to substitute derivatives U yy . In this way, obtain the equation: 

10a : 3 t / 4 + 156a; 2 x (1 + x U x ) U 2 + 36a; 2 U x U xx 

- 126a; 3 t / 4 - 126a; U 2 + 18a; U xx -18 U x - 252a; 2 C / 3 = 0 (17) 

This is a polynomial equation of degree 4 for U y , and all coefficients are functions 
of x only. By lemma [Gj this means U yy = 0, so U y = const =: c. 

For the second part of the proof, we insert this result into the p 4 -component 
of (llOal) . If we substitute U xx with the help of the Ernst equation, we find: 

6 U x (1 + x U x ) (1 + 2x U x ) = 0 

Hence, there are 3 cases: U x =0, U x = — ^ and U x = — We treat them 
separately: 

• If U x = 0 , use again the p 3 p 2 -component which reads 

10 a ; 3 c 6 = 0 , 


so c = 0 . 

• For f/j, = -j we have the same equation, so again c = 0. 

• The case U x = — y is slightly more involved. The p 3 P 2 -component reads 

c 2 (9 — 312a; 2 c 2 + 80x 4 c 4 ) 

8 a: 

Now, either c = 0 directly, or 9 — 312a; 2 c 2 + 80a; 4 c 4 . In the latter case, 
xc = const and hence c = 0 . 

□ 


Lemma 8 . Using Lewis-Papapetrou coordinates {x,y), assume the potential 
function U to he 


2 u 


= ku 


2 y + c + \J 4a ; 2 + 4 y 2 + Acy + c 2 


, withku,cG 1R. 


This provides a parametrization of flat space. 

Proof. Determine the function 7 from the secondary Ernst equations and find 


e 27 = 2 fc 7 y/ 4x 2 + 4 y 2 + 4cy + c 2 e 2U , fc 7 £ R. 

Then compute the Riemann curvature tensor for the metric 

g = e 2U (e -27 (dx 2 + dy 2 ) + x 2 dcf 2 ) - e~ 2U dt 2 

and find that it vanishes identically. Thus, the potential function U defines a 
flat metric, which of course is StAV. □ 
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Lemma 9. Let rkAl = 2 and ui = 0, but ^ 0. Assume a ^ 0. 

(ZTien t/iere is no non-trivial Killing tensor of valence 3. 

Proof. We assume there was such a Killing tensor. Then, by the necessary 
criterion (lemma [5]), A y — B x = 0. In addition, consider (llOal) . and the Ernst 
equations. Since we chose Lewis-Papapetrou coordinates with R = x, we have 

Uy ± 0. 

Consider (II Oal) in combination with the necessary criterion from lemma [5J 
plus the Ernst equations. The Ernst equations are to be invoked mainly in order 
to substitute (jp. We take derivatives w.r.t. x and y of (llOal) . Then, we have 
18 equations f (llOal) plus the necessary criterion A y — B x = 0 from lemma 0. 
Using the Ernst equations, we have only the following unknown functions: 

Uxxxxt bJxxxyi U X xx-> bJ xxyi Uxxi bJ xyi U X i U y , U, q. 

Use the ^-derivative of the pfp 2 -component to substitute U xxxy , and the x- 
derivative of the pfp|-component to substitute U XXX x dr terms of lower order 
derivatives. The quantity U xxy can be substituted for via the ^-derivative of 
the integrability criterion, but only if 

{l + 2xU x ){xUl-2,xUl + U x )^0. (18) 

In this case, we can proceed as follows: Substitute U xxx by the x-derivative of 
thep|-component, and use this component to substitute U xx . Finally, substitute 
U xy using the integrability condition. 

With all these substitutions at hand, we now have only equations in the 
unknowns U x , and U y left. For instance, the derivative w.r.t. y of the pf- 
component of (II Pall reads 

x Ul (1 + 2x U x ) (1 + x U x ) 2 (x U 2 X + U x + x U^f = 0. 

Therefore, either U x — P or U x = —^ or U x = or x Uf + U x + x U y = 
0. The three cases mentioned first are covered by lemma 0 and obviously in 
contradiction to the hypothesis U y ^ P. 

We are left with the forth case. Solve the equation and obtain 

U* = --U x (l + xU x ). (19) 

y X 

Then substitute this into the pf-component of (llPal) and obtain an expression 
for U xx , and from the integrability condition we obtain an expression for 
The other components of (llPal) are then satisfied trivially. At this point, it is a 
good idea to go back to the expressions for U y and U xy . Combining both, we 
find the equation 

l u ’ = - ixU Z 

1 One might want to check that both expressions are compatible, which is true. 
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which is an ODE for U y and can be solved in a straightforward way. The solution 
is 

1 

y V( Ax2 - h(y)) 

Use this to replace U 2 in (TTT71) : 


h(y) = -■ 


c(l + AxU x + Ax 2 U 2 ) 


U x (l + xU x ) 

Solve this for U x . There are two branches of possible solutions: 

JT _ 1 ~4x 2 - fi± sf Ax 2 /1 + fl 
x 2 x{Ax 2 + / 1 ) 

We can use the integrability criterion to find an explicit form for f\. First, 
obtain two differential equations: 


(fi) y 4\/7i = 0 


Up to the sign of the integration constant, both have the same solution 
fi = (2 y + c ) 2 = 4 y 2 + Acy + c 2 . 

Using again the equation for U 2 in terms of U x and integrating, one finds 

U = i ln(2y + c + \JAx 2 + 4 y 2 + Aye + c 2 ) + f 2 {x) 

with / 2 first being an unspecified integration ‘constant’. Checking if this solution 
is compatible with the expression for U x found above, ( f 2 ) x can take two possible 
values: Either ( f 2 ) x = -i, or 

— = (U) (x) = _^_ 

dx ( 2 y + c + a jAx 2 + ( 2 y + c) 2 ) \JAx 2 + ( 2 y + c ) 2 

Now, consider the Ernst equation U x + x U yy + x U xx = 0. For the first solution 
for (f 2 ) x , this implies x = 0, so this is no valid solution. Therefore, conclude 

U = i ln(2y + c + \J 4a: 2 + Ay 2 + Aye + c 2 ) — ln(a;) + c 2 

with an additional integration constant c 2 £ R. By lemma [SJ the metric is flat 
and therefore all Killing tensor fields are reducible. This concludes step (vi) in 
the list of section [ 2 ] 

To complete the proof, we still have to account for the case when (fTH]) is not 
satisfied. In this case, either U x = — ^ (this is covered by lemma 0 or 

U x (1 + x U x ) — 3xU 2 = 0 
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We solve for Uy\ 


j T 2 _ U x (1 + X U x ) 

y ~ 3x 


From the pf-component and the integrability criterion, we can also get another 
expression for Uy\ 


jj2 _ 3 U x (1 + x U x ) 


The only way to allow both solutions to be true is if U x = 0 or U x = — i. Both 
cases are covered by lemma [7) □ 

We have considered odd-parity third-degree integrals in a StAV space-time. 
Let us now summarize the results and merge them into one theorem: 

Proof of theo rem 1. If M is flat on a neighborhood, then it is totally reducible 
there Tho86). Thus assume that M is non-flat. 

First, consider only odd-parity integrals: 

Claim: Let M be non-flat with u> = 0, V V^ ^ 0, and V V u ^ 0. Let F be 
a third-degree integral of odd-parity in M. Then F is reducible by at least one 
degree. 

Proof of the claim. First, let us consider the case when there is an additional 
Killing vector field. As we have seen in lemma[3l this means that the odd-parity 
third-degree integral is reducible by the (linear) integral p y . So, the assertion is 
proven in this case. 

Second, if there is no additional Killing vector field, proposition 0 tells us 
(provided a/0) that there is no odd-parity third-degree integral. In the case 
a = 0, we have F = 0. Thus, the assertion is proven. □ 

Now, consider the even-parity contributions. The quadratic contributions 
( 2 ) ( 2 ) 

F) and Ff ' must obey the equation 


(T,if } } = 0 

as well as equations of the form 

{T,F^ k } + {V ab ,Fi 2) } = 0 , 


where a, &, k G {<f>, t}. 

These, however, are precisely equations msD for quadratic integrals with 
( 2 ) ( 2 ) 

leading term F ^ or Ff , respectively. This means that 

Pa ( F f) + F^ pi + f£ 1 p t Pt + F<$ p 2 t ), 

a = <f>,t, are quadratic integrals and therefore, the even-parity contributions to 
the degree 3 integral F are reducible by p^ and pt , respectively. Hence, also 
the entire integral F (consisting of the parts with degree from 3 down to 0) is 
reducible by one degree. □ 
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3.1 Zipoy-Voorhees 

Consider again the Zipoy-Voorhees class. We already considered third-degree 
odd-parity integrals in such space-times. Let us now consider even-parity com¬ 
ponents and assume w.l.o.g. 5^0. We use for the Hamiltonian H the repre¬ 
sentation 

H = Cl 1 p 2 x +n 2 pl + V < j > pl + VtPt (20) 

and denote the integral by 

F = a 0 p 2 x + ai p x p y + a 2 p 2 y + b 0 p\ + 6i p $ p t + b 2 p 2 . 

From each polynomial of degree 1 after split w.r.t. p x ,p y (cf. remarkQ]), we 
obtain integrability conditions for bo and b 2 . Automatically, b± = const is no 
longer of interest. 

Combining the Bertrand-Darboux relations and equations obtained from the 
degree 3 polynomial after splitting w.r.t. (p Xl p y ), we can solve for derivatives of 
do, ai and a 2 , and derive integrability conditions for them. From the integra¬ 
bility conditions, we can deduce that di = 0 and that (at least if S ^ 1) 

(y 2 - l)fl2 + ( x 2 - l)a 0 = 0. 

From the Bertrand-Darboux equations for b 0 , b 2 , we can now deduce d(do) in 
terms of do and solve the system of differential equations, obtaining 

d 0 = Cl (y 2 - x 2 y~ 52 (x + 1 f +s ~\x - l) 52 - 5 - 1 . 

Then, we can immediately compute a 2 : 

( x 2 — 1 \ S fx + l\ S x 2 — y 2 
a2 Cl \x 2 — y 2 ) \x — 1J y 2 — 1 

Finally, from the equations obtained from the degree 1 polynomial after split 
w.r.t. p x ,p yi we obtain the derivatives of 6 0 , b 2 , and by integration 

bo = -ci {y 2 - 1) {x 2 -1) + c 2 

62 = (jtt) + C3 

Comparing this result to the Hamiltonian shows that 

F = d H + c 2 p 2 ^ + c 3 p^ p t + c 4 p 2 . 

This means that every quadratic integral is reducible, provided 5 ^ 1 (in case 
S = 1, we obtain the Schwarzschild metric, which is integrable with the addi¬ 
tional integral in involution being given by a quadratic integral. This quadratic 
integral, however, is reducible by linear integrals that are not in involution). 
Together with theorem [l] this proves the assertion. □ 
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4 Conclusion 


In this paper we gave a proof for the reducibility of valence 3 Killing tensor fields 
in static and axially symmetric vacuum space-times (Weyl’s class). We saw that 
using prolongation-projection is an efficient way to decide on the existence of 
integrals in SAV metrics even if the metric is not given specifically. We plan to 
extend the result for degree 3 to the fully stationary SAV case. Though com¬ 
putationally more challenging, we are not aware of major conceptual problems 
with this. As for generalizations beyond the SAV context, the line of reasoning 
made here is in principle not specific to the SAV class of space-times, and an 
analogous approach might work for other classes of space-times, too. 
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